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Abstract 

The deformation kinetics for glassy polymers confined in microscopic domain at very low tem- 
perature regime was investigated using a transition-rate-state dependent model considering the 
shear thinning behavior which means, once material being subjected to high shear rates, the vis- 
cosity diminishes with increasing shear rate. The preliminary results show that there might be 
nearly frictionless fields for rate of deformation due to the almost vanishing shear stress in micro- 
pores at very low temperature regime subjected to some surface conditions : The relatively larger 
roughness (compared to the macroscopic domain) inside micropores and the slip. As the pore size 
decreases, the surface-to-volume ratio increases and therefore, surface roughness will greatly affect 
the deformation kinetics in micropores. By using the boundary perturbation method, we obtained 
a class of temperature and activation energy dependent fields for the deformation kinetics at low 
temperature regime with the presumed small wavy roughness distributed along the walls of an 
cylindrical micropore. The critical deformation kinetics of the glassy matter is dependent upon the 
temperature, activation energy, activation volume, orientation dependent and is proportional to the 
(referenced) shear rate, the slip length, the amplitude and the orientation of the wavy-roughness. 

PACS numbers: 83.60. St, 83.60.Rs,83.50.Lh 
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I. INTRODUCTION 



Glasses are amorphous materials of polymeric, metallic, inorganic or organic type. The 
plastic deformation of amorphous materials and glasses at low temperatures and high strain 
rates is known to be inhomogeneous and rate-dependent. In fact, the mechanical behavior 
of amorphous materials such as polymers [1-7] and bulk metallic glasses [8-10] continues to 
present great theoretical challenges. While dislocations have long been recognized as play- 
ing a central role in plasticity of crystalline systems, no counterpart is easily identifiable in 
disordered matter. In addition, yield and deformation kinetics [11-13] occur very far from 
equilibrium, where the state of the system may have a complex history dependence. 
In recent years, considerable effort was geared towards understanding how glasses respond 
to shear [14]. Phenomena such as shear thinning and 'rejuvenation' are common when shear 
deformation (rate) is imposed. At low temperatures they behave in a brittle elastic man- 
ner; at high temperatures, much above the glass transition the behavior is more (rubbery) 
(viscoelastic). There is a huge drop in modulus when the temperature is increased above 
the glass transition temperature, indicating a shift in behavior from (glassy) to (rubbery). 
Because of these peculiar mechanical properties of polymeric materials, the linear theory of 
viscoelasticity is unable to model closely the observed response and thus there is a need for 
a non-linear theory of viscoelasticity. 

Unlike crystals, glasses also age, meaning that their state depends on their history. When a 
glass falls out of equilibrium, it evolves over very long time scales. Motivated by the above 
issues and the interesting characteristics of deformation kinetics at very low temperature we 
shall study the deformation kinetics in microscopic domain at low enough temperature which 
is an interesting topic for applications in micro- and nanodomain [15-16] or the validation 
in using quantum mechanic formulations [9,17] where the nonlinear constitutive relations 
should be adopted. 

Meanwhile most of the classical solutions of contact problems, starting from the Hertzian 
case, rely on the assumption of nominally smooth geometries, which is reasonable at large 
enough scales. However, real surfaces are rough [18-19] at the micro- or even at the meso- 
scale, and the effect of multi-point contact is important for a series of phenomena involving 
friction and wear. The role of surface roughness has been extensively investigated, and op- 
posite conclusions have been reached so far. For instance, friction can increase when two 
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opposing surfaces are made smoother (this is the case of cold welding of highly polished met- 
als). On the other hand, friction increases with roughness when interlocking effects among 
the asperities come into play. This apparent contradiction is due to the effects of length 
scales, which appear to be of crucial importance in this phenomenon [15,20-21]. 
From the mechanical point of view, a contact problem involves the determination of the 
traction distributions transmitted from one surface to the other, in general involving normal 
pressures and, if friction is present, shear tractions, according to an appropriate set of equal- 
ities and inequalities governing the physics of the contact [22]. When there is friction at the 
contact interface, Coulomb friction behaviour is usually introduced to give the conditions 
necessary to determine the shear traction distribution. Any point in the contact area must 
be either in 'stick', or 'slip' condition, and the tangential tractions must behave accordingly. 
In this paper we shall consider the deformation kinetics of glassy polymer at very low tem- 
perature in micropores which have radius- or transverse-corrugations along the cross-section. 
The glassy polymer will be treated as a shear-shinning material. To consider the transport 
of this kind of glass (shear-thinning) polymer in microdomain, we adopt the verified model 
initiated by Cagle and Eyring [8] which was used to study the annealing of glass. To obtain 
the law of annealing of glass for explaining the too rapid annealing at the earliest time, 
because the relaxation at the beginning was steeper than could be explained by the bimolec- 
ular law, Cagle and Eyring [8] tried a hyperbolic sine law between the shear (strain) rate 
: £ and (large) shear stress : r and obtained the close agreement with experimental data. 
This model has sound physical foundation from the thermal activation process (Eyring [9] 
already considered a kind of (quantum) tunneling which relates to the matter rearranging by 
surmounting a potential energy barrier; cf., Fig. 1). With this model we can associate the 
(glassy) polymer with the momentum transfer between neighboring atomic clusters on the 
microscopic scale and reveals the atomic interaction in the relaxation of flow with (viscous) 
dissipation. 

The outline of this short paper is as follows. Section 2 describes the general mathemat- 
ical and physical formulations of the framework. In this Section, explicit derivations for 
the glassy deformation kinetics are introduced based on a microscopic model proposed by 
Eyring [9]. The boundary perturbation technique [18,23] will be implemented, too. In the 
third Section, we consider the very-low temperature limit of our derived solutions which 
are highly temperature as well as activation energy dependent at rather low temperature 
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regime. Relevant results and discussion are given therein. 

II. MATHEMATICAL AND PHYSICAL FORMULATIONS 

The beginnings of theoretical molecular mechanisms of deformation in amorphous poly- 
mers and glass are as old as the subject of atomic mechanisms of deformation and yield in 
metals. The first specific molecular mechanism of deformation for polymers and glass was 
published by Eyring [9] and later, Taylor [24] published the model of an edge dislocation to 
account for the plastic deformation in metals. However, whilst the theory of dislocations and 
crystal defects has become a major stream in the science of solid state, the corresponding 
effort applied to this problem in amorphous polymers must be considered rather small by 
comparison. 

The molecular theory of deformation kinetics came from a different stream of science than 
that of structure and motion of crystal defects (in particular dislocations). Its roots stretch 
to the developmental stages of theories of chemical reactions and thermodynamic descrip- 
tion of their temperature dependence, culminating in the key formulation by Arrhenius of 
the equation for reaction rates. By the beginning of this century the concept of activation 
entropy was included in the model, and it was considered that molecules go both in the 
forward direction (product state) and in the backward direction (reactant state). 
The development of statistical mechanics, and later quantum mechanics, led to the concept 
of the potential energy surface. This was a very important step in our modem understanding 
of atomic models of deformation. Eyring's contribution to this subject was the formal devel- 
opment of the transition state theory which provided the basis for deformation kinetics, as 
well as all other thermally activated processes, such as crystallisation, diffusion, polymeri- 
sation, etc. 

The motion of atoms is represented in the configuration space; on the potential surface the 
stable molecules are in the valleys, which are connected by a pass that leads through the 
saddle point. An atom at the saddle point is in the transition (activated) state. Under the 
action of an applied stress the forward velocity of a (plastic) flow unit is the net number 
of times it moves forward, multiplied by the distance it jumps. Eyring proposed a specific 
molecular model of the amorphous structure and a mechanism of deformation kinetics. With 
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reference to this idea, this mechanism results in a (shear) strain rate given by 

i V h k B T —AE . V h r 

where 

Vh — A2A3A, V m = A2A3A1, 

Ai is the perpendicular distance between two neighboring layers of molecules sliding past 
each other, A is the average distance between equilibrium positions in the direction of motion, 
A2 is the distance between neighboring molecules in this same direction (which may or may 
not equal A), A3 is the molecule to molecule distance in the plane normal to the direction 
of motion, and r is the local applied stress, AE is the activation energy, h is the Planck 
constant, k B is the Boltzmann constant, T is the temperature, Vh is the activation volume 
for the molecular event [9]. The deformation kinetics of the polymer chain is envisaged as 
the propagation of kinks in the molecules into available holes. In order for the motion of 
the kink to result in a plastic flow, it must be raised (energised) into the activated state 
and pass over the saddle point. This was the earliest molecular theory of yield behaviour in 
amorphous polymers, and Eyring presented a theoretical framework which formed the basis 
of many subsequent considerations. 
Solving Eqn. (1) for the force or r, one obtains: 

r^^sinh-Hl), (2) 

which in the limit of small (£/B) reduces to Newton's law for viscous deformation kinetics. 
We consider a steady deformation kinetics of the glassy polymer in a wavy-rough microtube 
of r Q (in mean-averaged outer radius) with the outer wall being a fixed wavy-rough surface 
: r = r Q + esin(k6) where e is the amplitude of the (wavy) roughness, and the wave number 
: k = 2n/L (L is the wave length). The schematic is illustrated in Fig. 2. Firstly, this 
material can be expressed as [9,18] £ = £ sinh(r/r ), where £ is the shear rate, r is the shear 
stress, and 

• 2k B TV h -AE 

^o = B= 7 T~T~ exp(— — , (3) 
h V m k B T 

is a function of temperature with the dimension of the shear rate, 

2k B T 

ro - — (4) 
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is the referenced (shear) stress, (for small shear stress r « r , the linear dashpot constitutive 
relation is recovered and r /£o represents the viscosity of the material). In fact, the force 
balance gives the shear stress at a radius r as r = —{r dp/dz)/2. dp/dz is the pressure 
gradient along the tube-axis or z-axis direction. 

Introducing the forcing parameter II = — (r o /2r )dp/dz then we have £ = £ sinh(Hr/r ). 
As the (shear) strain rate is 

(u is the rate of deformation (or velocity) in the longitudinal (z-) direct ion of the microtube), 
after integration, we obtain 

u = u s H — — ^ [cosh II — cosh ( — )], (6) 

11 Tq 

here, u s is the rate of deformation or velocity over the surfaces of the microtube, which is 
determined by the boundary condition. We noticed that a general boundary condition for 
material deformation kinetics over a solid surface was proposed (cf., e.g., [18]) as 

6u = L° a i(l-i)-^, (7) 

where 5u is the rate of deformation (or velocity) jump over the solid surface, L° s is a constant 
slip length and £ c is the critical shear rate at which the slip length diverges. The value of 
£ c is a function of the corrugation of interfacial energy. We remind the readers that this 
expression is based on the assumption of the shear rate over the solid surface being much 
smaller than the critical shear rate of £ c . £ c represents the maximum shear rate the material 
can sustain beyond which there is no additional momentum transfer between the wall and 
material-flow molecules. How generic this behavior is and whether there exists a comparable 
scaling for polymeric or amorphous materials remain open questions. 

At small pressure gradient, the shear-thinning matter behaves like a Newtonian flow, while at 
high pressure gradient, the shear-thinning matter flows in a plug- flow type. Such a behavior 
is due to the shear thinning of the material, i.e., the higher the shear rate is, the smaller is 
the (plastic) flow resistance [7]. On the microscale, this shear-thinning matter can bridge 
the Newtonian deformation kinetics to that of the pluglike type and offers us a mechanistic 
model to study the deformation kinetics in micro- and even nanodomain using the technique 
of continuum mechanics. 
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With the boundary condition from (cf., e.g., [18]), we shall derive the rate of deformation 
(or velocity) field or deformation kinetics along the wavy-rough microtube below using the 
boundary perturbation technique (cf. [23]) and dimensionless analysis. We firstly select the 
hydrodynamical diameter L r to be the characteristic length scale and set 

r' = r/L r , R = r /L r , e' = e/L r . (8) 

After this, for simplicity, we drop all the primes. It means, now, r, R Q , and e become 
dimensionless. The wall is prescribed as r = R Q +e sm(k9), and the presumed fully-developed 
plastic flow is along the z-direction (microtube-axis direction). Along the confined (wavy) 
boundaries, we have the strain rate 

£ = ^rfn^ on surface' (9) 
where, n means the normal. Let the rate of deformation u be expanded in e : 

u = uq + eui + e 2 u 2 H , (10) 

and on the boundary, we expand u(r + edr, 9{= 6q)) into 

dv 2 

u(r,9)\ {ro+£dri e o) = u(r ,9) + e[dru r (r , 9)] + e 2 [—u rr (r , 9)] H = 

{u slip + ^ coshffi|^+- n W ro = R o] (11) 
where the subscript means the partial differentiation (say, u r = du/dr) and 



siiplon surface on surface 

(12) 



u slWo = ^o[sinhn(i - k^Eyv^ (13) 



Now, on the outer wall (cf., e.g., [23]), the (shear) strain rate 
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• du „ V(r — R — esm(k9)) r , ~k , ., , 

k 

^cos(k9)u e \ {Ro+tdrfi) ] =u 0r \ Ro + e[u lr \ Ro + u 0rr \ Ro sm(k9)- 



r 

k lk 2 
r -2«o 9 k C0S (^)] + e 2 [-2^ 2 



-oMoelflo cos(fc0)] + e 2 [--— cos 2 (k9)u 0r \ Ro +u 2r \ Ro +u lrr \ Ro sin(k9) + 



1 k 

-u 0rrr \ Ro sin 2 (k9) - -cos(k9)(u le \ Ro + u 0er \ Ro sm(k9))} + 0(e 3 ). (14) 



Considering L° ~ R Q ^> e case, we presume sinhll <C £ c /£o so that we can approximately 
replace [1 - (f„ sinh n)/4]~ 1/2 by [1 + £ sinh n/(2£ c )]. With equations (6), (7), (9), (10), 
(11) and (14), using the definition of the (shear) strain rate £, we can derive the rate of 
deformation (or velocity) field up to the second order. The key point is to firstly obtain the 
slip rate of deformation (or velocity) along the wavy boundaries or surfaces. 
After lengthy mathematical manipulations and using (1 — i/ic)" 1 ^ 2 ~ 1 +£/(2£ c ), 

m = -^^[cosh(^) - coshlTj + u sliPo , (15) 



Mi = ^ o sm(k0) sinhll + u sUpi , (16) 



we have 



U slip = m-U 0r (l - ^)}\r=R + 6[-«/(l - ^)}\r= Ro + ^ ~ U sc (l - ^)]\ r=Ro } = 



where 



Usiipo + eu sHpi + e u sHp2 + 0(e) (17) 



u 0r = -^osinh(^r), (18) 

u 0rr = -io^-cosh(^-r), (19) 

• n 2 n 

u 0rrr = -£ —smh( — r), (20) 



k n • n 

u f = u lr + u 0rr sin(k9) - — cos(k6)u 0e = -— £ cosh(— r) sin(A;0), (21) 



and 



u *° = cos2 (^) M o r + \u 0rrr sm 2 (ke) = \Ut^ cos 2 (k6) - ^ sm 2 (k6)] sinh(-^r). (22) 

Thus, at r = R Q , up to the second order, 

u sHp = u s = L% sinh n(l + ^) + e^o sm(k6) [sinh IT + ^-L° s cosh n (1 + K Q )\ + e 2 L° s ^{[ 

2 K I 

TT rnsh TT k 2 TT 2 TT 2 P 

p r0 sin 2 (A;0) - ^- cos 2 (A;0) + ^ sin 2 (^)] sinh n(l + K Q ) + cosh 2 n sm 2 {k6)}, 
R L S R R Q R Q £ c 

(23) 

where 

K = l + (i smhU)/i (24) 



From the rate of deformation (or velocity) fields (up to the second order), we can integrate 
them with respect to the cross-section to get the volume (plastic) flow rate (Q, also up to 
the second order here). 

r6p rR +(.s\ri(kd) 

Q = J Q J M ( r ' e ) rdrde = Q smooth + e Qp + e 2 Q P 2 • ( 25 ) 

III. RESULTS AND DISCUSSION 

We firstly check the roughness effect upon the shearing characteristics because there are 
no available experimental data and numerical simulations for the same geometric configu- 
ration (microscopic tubes with wavy corrugations in transverse direction). With a series of 
forcings (due to imposed pressure gradients) : II = R o (—dp/dz)/(2r ), we can determine 
the enhanced shear rates (d£/dt) due to forcings. From equation (5), we have (up to the 
first order) 

^ = ^[sinhn + esin(M)^coshn]. (26) 
dt dt rl 

The calculated results are demonstrated in Figs. 3 and 4. The parameters are fixed below 
(the orientation effect : sm(k9) is fixed here). r Q (the mean outer radius) is selected as the 
same as the slip length L° s = 100 nm. The amplitude of wavy roughness is e = 0.04, 0.07, 0.1, 
the Boltzmann constant (ks) is 1.38 x 10 -23 Joule/°K, and the Planck constant (h) is 
6.626 x 10~ 34 Joule • s. 

In each panel, the inner curve is the relevant boundary of the tube or the geometric part 
of the presentation. The distance between the inner and corresponding outer curves is the 
calculated physical shear rate : £. We can observe once the temperature (T) changes a little 
from 285 ° to 295 °, the enhancement of £ becomes at least three orders of magnitude (for 
II = 1, the activation energy : 3 x 10~ 23 Joules). Even at very low temperature Fig. 4 
gives very large strain rates which are required to to obtain the necessary strain for plastic 
deformation. Thus, the constitutive relations is highly nonlinear at rather low temperature 
regime [9]. It is worth pointing out that the Eyring model requires the interaction between 
atoms in the direction perpendicular to the shearing direction for the momentum transfer. 
This might explain why our result is orientation dependent. The effect of wavy-roughness 
will be significant once the forcing (II) is rather large (the maximum is of the order of mag- 
nitude of e[ntanh(n)/ J Rj). 
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To be specific, we can illustrate the shear rate (£) with respect to the temperature (T) once 
we calculate £ as the latter is temperature dependent (but presumed roughness independent 
here) which could be traced from equation (1). This is shown in Fig. 5. 
Note that, based on the rate-state Eyring model (of stress-biased thermal activation), struc- 
tural rearrangement is associated with a single energy barrier (height) E that is lowered or 
raised linearly by a (shear) yield stress r. If the transition rate is proportional to the plastic 
(shear) strain rate (with a constant ratio : Co; £ = CoR t , Rt is the transition rate in the 
direction aided by stress), we have r = E /V* + (k B T/V*) ln(^/C ^o) or 

"£^>*&- 

where V* = Vh is a constant called the activation volume, ks is the Boltzmann constant, T 
is the temperature, u is an attempt frequency or transition rate [8,25], and f is the stress 
rate. Normally, the value of V* is associated with a typical volume required for a molecular 
shear rearrangement. Thus, if there is a rather-small (plastic) flow (of the glass) at low 
temperature environment then it could be related to a barrier-overcoming or tunneling for 
shear-thinning matter along the wavy-roughness (geometric valley and peak served as atomic 
potential surfaces) in cylindrical micropores when the wavy-roughness is present. Once the 
geometry-tuned potentials (energy) overcome this barrier, then the tunneling (spontaneous 
transport) inside wavy-rough cylindrical micropores occurs. 

To examine the behavior of the shear rate at low temperature regime, we calculate £ and 
£ (CqVq = 5 x 10 10 s~ x ) with respect to the temperature T and show the results in Fig. 5. 
For a selected activation energy : 5 x 10~ 22 Joule or ~ 10~ 3 eV (a little bit smaller than 
the binding energy of 3 He), we can find a sharp decrease of shear rates around T ~ 0.01°K. 
Below this temperature, there might be nearly frictionless transport of glassy matter. Note 
also that, according to Cagle and Eyring [8], V* = 3V5£/2 for certain material during an 
activation event, where V is the deformation volume, 5£ is the increment of shear strain. 
If we select a (fixed) temperature, say, T = 0.1°K, then from the expression of r , we can 
obtain the shear stress r corresponding to above forcings (n) : 

r = r sinh _1 [sinh(n) + esin(fc0)-5- cosh(n)]. (28) 

-Re- 
There is no doubt that the orientation effect (9) is also present for deformation kinetics of 
polymeric matter. For illustration (shown in Fig. 6), we only consider the maximum case 
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: | sm(k9)\ = 1. The trend of enhancement due to II (pressure-forcing) and e (roughness) 
is similar to those presented in Figs. 3 and 4. We remind the readers that, due to the 
appearance of r , we fix the temperature to be the same and the activation volume : 1CT 25 
m . 

In fact, as shown in Fig. 6, the calculated (shear) stress (which is directly linked to the 
resistance of the glassy matter) also shows a sudden decrease around T ~ 1.1°K especially 
for the case of (Co^o = 2 x 10 9 s" 1 ). Here, the activation volume (V* or Vh) is selected as 
0.2 nm 3 [25]. Thus, the nearly frictionless transport of the glassy fluid at low temperature 
environment (relevant to the supersolidity, cf. [26]) could be related to a barrier-overcoming 
or tunneling for shear-thinning matter along the wavy-roughness (geometric valley and peak 
served as atomic potential surfaces) in cylindrical micropores when the wavy-roughness 
is present. Once the geometry-tuned potentials (energy) overcome this barrier, then the 
tunneling (almost frictionless transport) inside wavy-rough cylindrical micropores occurs. 
We also noticed that, as described in [9], mechanical loading lowers energy barriers, thus 
facilitating progress over the barrier by random thermal fluctuations. The simplified Eyring 
model approximates the loading dependence of the barrier height as linear. This Eyring 
model, with this linear barrier height dependence on load, has been used over a large fraction 
of the last century to describe the response of a wide range of systems and underlies modern 
approaches to sheared glasses. The linear dependence will always correctly describe small 
changes in the barrier height, since it is simply the first term in the Taylor expansion of 
the barrier height as a function of load. It is thus appropriate when the barrier height 
changes only slightly before the system escapes the local energy minimum. This situation 
occurs at higher temperatures; for example, Newtonian deformation kinetics is obtained 
in the Eyring model in the limit where the system experiences only small changes in the 
barrier height before thermally escaping the energy minimum. As the temperature decreases, 
larger changes in the barrier height occur before the system escapes the energy minimum 
(giving rise to, for example, non- Newtonian deformation kinetics). In this regime, the linear 
dependence is not necessarily appropriate, and can lead to inaccurate modelling. This 
explains why we should adopt the hyperbolic sine law [9] to treat the glassy matter. 
To be specific, our results are rather sensitive to the temperature (T) and the activation 
energy. Fig. 7 shows especially the temperature dependence of the forcing parameter (IT) 
if dp/dz is prescribed (say, around 6 x 10 10 Pa/m) and the activation volume is 0.2 nm 3 
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(r = 100 nm). We can observe that once T increases II decreases. £ calculated using 
prescribed II and using directly T also differs. Finally, we present the calculated maximum 
velocity (unit : m/s) with respect to the temperature in Fig. 8. Geometric parameters : 
r Q and the activation volume are the same as those in Fig. 7 and the roughness amplitude 
e = 0.02, 0.05i? o . We consider the effect of the activation energy : 9.0 x 10 -23 , 1 x 10~ 22 
and 2 x 10~ 22 Joule. Around T ~ 0.35°K, the maximum velocity (of the glassy matter) 
either keeps decreasing as the temperature increases for larger activation energy or instead 
increases as the temperature increases for smaller activation energy! The results presented 
in Fig. 8 might be related to the microscopic origin for physical aging or effects of thermal 
history [27] and indeterminate solutions discussed in [28]. The latter observation might be 
related to the argues raised in [29] for the annealing process of solid helium at similar low 
temperature environment if we treat the solid helium to be glassy at low temperature regime. 



IV. CONCLUSIONS 



To give a brief summary, we analytically obtain a class of temperature as well as 
activation energy dependent fields of the rate of deformation for glassy polymer in 
microscopically confined wavy-rough domain at very low temperature regime. The effects 
of wavy corrugation upon the confined deformation kinetics at very-low temperature are 
clearly illustrated. It is found that there exist almost frictionless plastic flow fields for the 
rate of deformation of glassy polymer inside cylindrical micropores at very low temperature 
once the micropore surface is wavy-rough and the activation energy is prescribed. Once 
the temperature, activation volume, and geometry are fixed, the increase of activation 
energy instead reduces significantly the (maximum) rate of deformation of the glassy 
matter. The critical rate of deformation is proportional to the (referenced) shear rate, 
the slip length, the orientation and the amplitude of the wavy-roughness as illustrated above. 
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Fig. 1 Increasing strain causes a local energy minimum to flatten until it disappears 
(removing of energy barrier or quantum-like tunneling) . The structural contribution to 

the shear stress is shear thinning. 
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Fig. 2. Schematic (plot) of a micropore, e is the amplitude of small wavy-roughness. 
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Fig. 3. Comparison of the shear rate (£) of polymeric matter in smooth microtubes and 
wavy-rough microtubes for k = 10, e = 0.0,0.04,0.07,0.1 with r o =100 nm. Cc/Co = 10 and 
= r Q . k is the wave number and e is the amplitude of the wavy-roughness. T is the 
temperature. The solid-line length represents the scale of ( = 500s -1 . 
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Fig. 4. Comparison of the shear rate (£) of polymeric matter in smooth microtubes and 
wavy-rough microtubes for k = 10, e = 0.0,0.04,0.07,0.1 with r o =100 nm. Cc/Co = 10 and 

= r Q . k is the wave number and e is the amplitude of the wavy-roughness. T is the 
temperature. The solid-line length represents the scale of £ = 500s -1 . As the temperature 

increases a little, £ increases significantly. 
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Fig. 5 Comparison of calculated shear (strain) rates using an activation energy 5 x 10 22 
Joule or ~ 1CT 3 eV. There is a sharp decrease of shear rate around T ~ 1.4°K. 

C u = 5 x lO 10 *- 1 (u ~ 12*). 
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Fig. 6 Comparison of calculated (shear) stresses using an activation energy 5 x 10 22 
Joule or ~ 1(T 3 eV. There is a sharp decrease of shear stress around T ~ 1.1°K for 
Co^o = 2 x 10 9 s _1 . Below 1.1°K, the transport of polymeric matter is nearly frictionless. 
uq is an attempt frequency or transition rate 25 . 
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Fig. 7 Calculated forcing parameters (x/-R ) w.r.t. the temperature (T). 
ftsoM is the compressibility, ^4 is the area, and V2 is the selected volume. 
Forcing (%) decreases as the temperature (T) increases. 
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Fig. 8 Comparison of calculated (maximum) velocity (unit : m/s) using two activation 
energies 1.5 x 10~ 22 and 2 x 10 -22 Joule. Around T ~ 0.35°K, the monotonic trend of 
velocity (or deformation rate) bifurcates as the temperature increases. 
r a = 100 nm and e = 0.02, 0.05i? o . 
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